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Dwupasmowy model metalu w emisji polowej elektronów 
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Two-Band Model in Field Emission of Electrons from Metals 


WSTĘP 


Obliczenia struktury pasmowej metali przejściowych stanowią 
bardzo skomplikowany problem. Struktura pasmowa tych metali odbiega 
bowiem zdecydowanie od struktur, jakie posiadają metale proste, Wię= 
kszość metod obliczeń struktury pasmowej ciał stałych można podzie- 
ué na dwie grupy. Jedna z nich obejmuje metody polegające na diago- 
nalizowaniu jednocząstkowego hamiltonianu wyprowadzonego z pierw 
szych zasad. Sa to np.: metoda funkcji Greena (KKR), metoda uzupeł- 
nionych fal płaskich (APW). Druga grupa zawiera metody polegające 
na dopasowywaniu hamiitonianu w taki sposób, aby otrzymane wartoś- 
ci cneraii dobrze odtwarzały realne struktury, Do nich można zaliczyć 
wszelkiego rodzaju półempiryczne metody, bazujące na modelowych 
pscudopotencjatach, oraz metody interpolacyjne, proponowane np. 
przez Slatera i Koste ra [1] czy Mueller af2]. 
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Jednym z prostszych podejść do problemu metali o skompliko= 
wanej strukturze energetycznej jest założenie, że elektrony "d" i ele- 
ktrony "s" w tych metalach są "różnymi" cząstkami. Oddziaływanie 
między nimi prowadzi do przechodzenia jednego typu cząstek w dru- 
gie - rozpraszanie s-d, Zdejmuje ono zarazem degenerację wartości 
własnych energii pasm "s" i pasm "d" odpowiedniej symetrii w punktach, 
gdzie pasma te przecięłyby się w nieobecności takiego oddziaływania, 
W metalach przejściowych poziom Fermiego znajduje się właśnie w 
obszarze, w którym zachodzi dosyć mocne mieszanie się stanów typu 
"s" i typu "d", Prowadzi to do znacznego odejścia pasm ehetgetycx- 
mych od swojej, nie zaburzonej oddziaływaniem postaci, Wszystkie 
dotychczasowe rachunki, dotyczące emisji polowej elektronów z meta- 
li, nie uwzględniały wpływu istniejących zhybrydyzowanych pasm na 
charakterystyki emisyjne. Ponieważ, jak już wapomnieliśmy, stany le- 
zace w pobliżu energii Fermiego są szczególnie zdeformowane przez 
oddziaływanie sed, należy spodziewać się, że zaniedbanie tych efek- 
tów może zniekształcać niektóre charakterystyki emisyjne. 

W niniejszej pracy podamy opis emisji polowej elektronów z mo= 
delowo zhybrydyzowanych pasm energetycznych, Pasma energetyczne 
emitera otrzymamy, przyjmując modelowy Sinusoidalny potencjał krys- 
taliczny. Taki wybór potencjału uwarunkowany jest żądaniem, aby 
potencjał w miarę dobrze odtwarzał realny potencjał krystaliczny oraz 
prowadził do analitycznego rozwiązania równania Schrbdingera. W 
rozdziale drugim podamy rozwiązanie równania Schrédingera w przy 
padku mocelowej hybrydyzacji pasm energetycznych. W rozdziale 
trzecim rozpatrzymy rozkłady energetyczne emitowanych elektronów z 


pasm zhybrydyzowanych I porównamy je z rozkładami otrzymywanymi 
w modelach swobodnoelektronowych, 


ENERGETYCZNY MODEL EMITERA 


W rozdziale tym określimy strukturę energetyczną ! funkcje falo- 
we elektronu w krysztale przy modelowo wprowadzonej hybrydyzacji 
nałożonych na siebie pasm, Zakładamy, że w krysztale istnieją dwie 
grupy elektronów. Jedna z nich opisana jest modelowym potencjałem 
sinusoidalnym o takiej amplitudzie, aby dozwolone drugie pasmo ener 
getyczne było możliwie szerokie. Pasmo to będzie odgrywało rolę 
swobodnoelektronowego pasma "s". Druga grupa elektronów opisana 


jest także modelowym sinusoidalnym potencjałem krystalicznym, ale 
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o amplitudzie odpowiednio większej, Prowadzi to do pasm wąskich i w 
sposób przybliżony pozwala modelować pasmo mocnego wiązania [3]. 
Amplitude potencjału krystalicznego obu grup elektronów dobieramy tak, 


aby wąskie pasmo "d" znajdowało się wewnątrz szerokiego pasma "s", 
Dla każdej z grup elektronów wprowadzamy potencjał postaci: 


V(z)=2W, (1+cos 2) , z<0 (1A) 


gdzie: GW, - amplituda potencjału, 
g - stała sieci kryształu, 


Rozpatrując w dalszej części pracy emisję elektronów, przyj- 
miemy dla z > O potencjał w postaci: 


V(z)=-eFz -E, + (1B) 


gdzie: F ~- natężenie pola elektrycznego, 
Ept- energia Fermiego i praca wyjścia elektronu z metalu. 


Ryc. 1. Model zhybrydyzowanych pasm ener- 
getycznych; szerokość pasm przed hybrydy= : e 
zacją 11 eV i 1 eV, szerokość przerwy ener- 02 z zal 

getycznej 0,15 eV kla) (K* ) — 


Funkcję falową elektronu dła wnętrza emitera zapiszemy w 


postaci: 
dż 
oe eee (2) 
-+ - > 
gdzie: ka I, ~ składowe wektorów k p równoległe do powierz- 


, 
chni emitera. 


Taka postać funkcji falowej uzasadniona jest jednowymiarowym charak 
terem bariery przypowierzchniowej, Ponieważ w płaszczyźnie prosto- 
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padłej do kierunku emisji nie działają żadne siły, składowa wektora 
falowego elektronu, równoległa do powierzchni metalu, zostaje w trak- 
cie tunelowania zachowana, 


Równanie Schridingera transformuje się więc do postaci: 


- A — +(V(z)-E,) (z)=0 (3) 


2,2 
gdzie: energiaW= A - energia "normalna", tzn. ener- 
| gia elektronu związana ze składową pędu w kierunku emisji. 

Przyjmując za funkcję falową elektronu, funkcję wyrażoną wzo- 


rem (2) otrzymujemy: 


a+ (a +24c0825) 9(8) =0 (4) 
gdzie: 3 = am fh Wis 


Rozwiązanie powyższego równania zapiszemy w postaci; 


o P 
PE Dabare złe) 
ej (5) 
„pźłz = -LEz(2r 
BE > Cz. © ‘ e 


r. 


gdzie: 1 = 0,1 odnosi sie do parzystych i nieparzystych dozwolonych 
obszarów energii elektronu w krysztale, 
A, B - stałe całkowania, 


p Cop} 7 WSpółczynniki podlegające wyznaczeniu, 
Zauważmy, że postać funkcji (5) ma kształt ogólnej reprezen- 


tacji funkcji falowej w krysztale: 


9(z) m eug (z) Be "u, z (2) . ( 6 ) 
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gdzie: N - numer pasma, 
Ung (Z) - periodyczna część funkcji falowej. 
We wzorze (5) pierwszy (drugi) wyraz odpowiada fali elek- 


tronowej padającej (odbitej) na barierę powierzchniową od strony me- 
talu. 


Rozwiązanie równania Schrédingera w obszarze z>O ma pos- 


= VJ 
dia  z)z,, z," (E-*P-E,)/eF, v- (2996 (2-2), 


f(z) “A e wyk Da (4 v») 


(7A) 
+1,(55)87%] 
da  Z<z,, YacY, 


piar r gY) 2 


-1,(¢y™")e*] m 


gdzie: Jy (x) „ly (x) - funkcje Bessela; pierwszego rodzaju, zmo- 
dyflkowana 
Wprowadémy dalej modelowy potencjał oddziaływania elektro- 
nów sed Voa (k) « Niech E, (k) , E (k) oznaczają wartości 


własne hamiltonianu nie oddziaływających elektronów. Zhybrydyzowane 
pasma otrzymamy poprzez diagnonalizację modelowego hamiltonianu: 


He E(k) Y) 
W EK 


Va Ü) p Vs Q 


(8) 
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Wartościami własnymi powyższego hamiltonianu są: 


+ s 2) 2 
EOF e Ge EA- | (s) 


gdzie +(-) oznacza górne (dolne) zhybrydyzowane pasmo energetycz- 
ne. Niech f lk) i g_(k) będą najni falowymi nie oddziaływających 
pasm "s" i "da", Funkcje falowe g (f) pasm zhybrydyzowanych wyra- 
żają się poprzez kombinacje funkcji (0 i Sy(K) w sposób następują- 
cy: 


4 Hg @ gre ga) (20) 


gdzie: aż R), a (k) są współczynnikami hybrydyzacyjnymi fun- 
kcji pasm nie odziaływających i np, dla pasma górnego Spełniają zwia- 
zki (podobne związki spełniają i dla pasma dolnego): 


a, (X) Jaz *v,„(WY/(EG)-E,(w)) (11) 


zakładając dalej, że funkcje pasm nie oddziaływających są 


unormowane i ortogonalne, otrzymujemy (dla pasma górnego): 


0,(R)-[(E,(®)-€())/(E, ()-€,(8)-2€ (5) |” (12) 
ARKE A-E EE 0 -2E() ] 


Obecnie przejdziemy do konstrukcji funkcji falowych. Jak już 
wspomnieliśmy, pasmo "s" modelujemy poprzez użycie potencjału o 
dosyć małej amplitudzie ~ 1 eV. Otrzymamy pasma energetyczne 
o szerokościach odpowiednio: pierwsze pasmo = 3,5 eV i drugie ~ 
11,5 eV. Jako funkcje typu "s* weźmiemy funkcje falowe odpowiadające 
drugienu pasmu. Pasmo to, zarówno szerokością, jak i charakterem 
funkcji falowych odpowiada większości pasm swobodnoelektronowych 
spotykanych w realnych metalach. Pasmo "d" modelujemy, przyjmując 
potencjał krystaliczny także w postaci sinusoidalnej, ale o amplitudzie 
odpowiednio dużej ~ 19 eV. Prowadzi to do powstania pierwszego 
pasma energetycznego o szerokości ~ 1,8 eV, znajdującego sie wew- 
nątrz szerokiego pasma "s". Jeszcze raz przypominamy, że z założe- 
nia przyjmujemy, iż obie grupy elektronów sa niezależne. Każda gru- 
pa clektronów "widzi" tylko "swój" potencjał, funkcje falowe zaś z 
założenia są ortogonalne. Potencjał oddziaływania V alk) traktujemy 
również modelowo. Jego kształt przyjmiemy w postaci podobnej do za- 


leżności dyspersyjnej pasm nie oddziaływających. 
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Tak więc, Veal’) przyjmiemy w postaci takiej zależności 
energii od wektora falowego, jak w przypadku potencjału sinusoidalne- 
go o pewnej amplitudzie. Amplitude tę dobieramy tak, aby pasma po 
zhybrydyzowaniu nie odbiegały zbytnio na granicach strefy Brillouina 
od swoich położeń pierwotnych (małe przesunięcia hybrydyzacyjne). 
Kształt pasm przed i po zmieszaniu przedstawia ryc. 2. 


Ryc. 2. Model zhybry tyzowa- 
mych pasm energetycznych; -A “*, o ka % 
opis w tekście — 


Funkcja £00, odpowiadająca drugiemu pasmu, ma postać w 
schemacie strefy rozwiniętej (przyp. jednowym.): 


= fede, irz 
pikire wL 


a w schemacie strefy zredukowanej: 


pt ez = s Bfr) 
AGE 2 ać + (14) 


(13) 


j 2 s 
gdzie: B, 7478 r 
Funkcja gal). reprezentująca pasmo wąskie "d", wyraża się wzorem: 


d se = 
TOR Fata. n (35) 


Rozpatrzmy zhybrydyzowane dolne pasmo (ryc. 2). Niech Kin 
oznacza wektor falowy, odpowiadający maksimum tego pasma, Oznacza 
my obszar w przestrzeni wektora falowego od O do km numerem pier- 
wszym, od O do — km numerem drugim, od = km do — z numerem trze- 
cim Í od + km do + — numerem czwartym. Zauważmy, że stany 
elektronowe opisane wektorem falowym z obszarów pierwszego i trze- 
ciego przedstawiają fale elektronowe, poruszające się na prawo, Wy- 
nika to ze znaku pochodnej fh" dE/dk= Vz 20% Natomiast ele- 
ktrony o wektorach falowych z obszarów drugiego i czwartego poru- 


szalą się w kierunku na lewo. 
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ROZKŁADY ENERGETYCZNE EMITOWANYCH ELEKTRONÓW 
PORÓWNANIE Z MODELEM ELEKTRONÓW SWOBODNYCH 


Rozpatrzmy emisję z metalu o strukturze energetycznej przed- 
stawionej na ryc. 1. Funkcję falową elektronu poruszającego się w 
metalu zapiszemy w postaci: 


METAL 1-0 =», 
4  (z):A%(z)+BY(z) (16) 
gdzie strzałki oznaczają kierunek rozprzestrzeniania sie fal elektronc= 


wych. W przypadku, gdy wektor falowy elektronu znajduje się w ob- 
szarze pierwszym (ryc. 2), to (pasmo dolne) 


¥ (z) = a eT) ct Ce Fraz, +0; (B)e ipft Pe 


Pa-eo 


(17) 


Gdy wektor falowy elektronu znajduje się w obszarze trzecim, otrzy- 
mujemy: 


+ (z) = a, ( pe tis p e Pda, {B) Pale Bod ptre (18) 


re-a 


Fale elektronowe odbite od bariery otrzymujemy poprzez sprzężenie 
powyższych wyrażeń. 

Współczynnik przejścia elektronów z metalu przez powlerz= 
chniową barierę potencjału obliczymy, biorąc stosunek prądu elektrono- 
wego przechodzącego przez barierę powierzchniową i padającego na 
barierę od strony wnętrza emitera, W przypadku, gdy wektor falowy 
padającego elektronu znajduje się w obszarze pierwszym (ryc. 2), 


otrzymujemy: _ srj 
47h sry pad 


(19) 


220 
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4, * Ee yc „(Pr2r) , 


an pice Ż_ cl (B+2r) ys > Chul- zz), 
[c/a]- 04,0, /| Gaye Bay | | 


onone) Chato CS, 


izao ) 


+ {af 6. przygo gea), 


gdzie: © (Z) wyraża się wzorem (7). 


W przypadku, gdy tunelujący elektron posiada wektor falowy z obsza 
ru drugiego (ryc. 2), we wzorze na współczynnik D zachodzą zmiany: 


Gy -— - G 

G; ——> - Go 

A Ap = A 
21 21 


Rozkład energetyczny wyliczymy z ogólnie stosowanego wzoru 


[3]: 


TED(E)* “Ra exp((E- wey | DE „JdE, (20) 


Obliczenia numeryczne przeprowadzono dla następujących para- 
metrów: ‘P = 4,5 eV, T = 76° K, szerokość pasma "s"~ 11 eV, sze- 
rokość pasma "d" -—1 eV, E- = 8,2 eV (względem dna pasma zone 
kiego przed hybrydyzacją). Wyniki obliczeń dla pól F = 2 ° 107 V/cm 
i6» 107 Vicm przedstawione są na ryc, 3. Dla porównania przedsta- 
wiano na niej również sumę rozkładów energetycznych z pasma sze- 
rokiego "s" | wąskiego "d". Pozkłady z pasm zhybrydyzowanych od- 
zwierciedlają dość ostro charakterystyczne punkty (A,B,C,D) struktury 


pasmowej (ryc.3). Punkt A (pojawienie się nowego pasma) daje zauwa- 


| „z RO TEE na i 


Ryc. 3. Rozkłady energetyczne 
emitowanych elektronów z met 
- - . lu o modelowo zhybrydyzowanej 
—E-€, lew) —— E-E, tev) strukturze pasmowęj 

żalny wkład do rozkładu, praktycznie tylko dla pola 6 «10° Vian. 
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Natomiast punkt B zaznacza się na obu krzywych dosyć wyraźnie, 
Związane tosjest ze zmniejszaniem się współczynnika przejścia w 
miarę zbliżania się do maksimum dolnego pasma, Krzywa rozkładu 
energetycznego elektronów pochodzących z górnego pasma również 
odzwierciedla punkty C i D struktury pasmowej, 

Rozkłady, począwszy od dna górnego pasma, Szybko rosną. 
W punkcie D daje się zauważyć pewne zaniżenie emisji, Związane 
jest to z prawym brzegiem górnego pasma. 

Dla wyższych energii rozkład nie wykazuje już żadnych 
osobliwości, Rozkład energetyczny elektronów emitowanych z metalu 
o strukturze przedstawionej na ryc, 1, Ale bez hybrydyzacji, przedsta 
wiałby sie jako suma rozkładów z poszczególnych pasm, natomiast 
w obszarze energetycznym, odpowiadającym przckrywanlu się tych 
pasm, charakteryzowałby się pewrym wzmożenien. emisji, Dodanie do 
tego modelu oddziaływania powodującego hybrydyzację pasm energe 
tycznych daje w wyniku przerwę energetyczną. Prowadzi to w prostej 
linii do zaniku prądu emisyjnego w tym obszarze energii, W doświad- 
czeniu, przy emisji zarówno z wolframu, jak i molibdenu, nie obserwu- 
jemy tak gwałtownych zmian w rozkładach energetycznych emitowanych 
elektronów. Wynika to stąd, że w strukturze energetycznej tych mete 
i, w obszarze energii badanych za pomocą emisji polowej, nie ma 
absolutnych przerw energetycznych. Ne przykład w wolframie, w 
strukturze energetycznej wzdłuż xicrunku 4 (w okolicy poziomu 
kermiego), istniejąca przerwa między pasmami o symetrii A Jes! prec- 
kryta dwoma pasmami: ô, i p. Gdybyśmy obliczyli rozkład energe- “ 
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tyczny elektronów emitowanych z metalu o takiej strukturze pasmowej 
(z założenia izotropowej) - tzn. wolfram, kierunek Q ~ to nie otrzy- 
malibyśmy tak ostrej przerwy jak na ryc.3. Pasma przekrywające przer- 
wę także dałyby bowiem swój wkład do emisji i efekt przerwy hybrydy= 
zacyjnej byłby częściowo ukryty.Niemniej jednak,rozkład ten odzwier- 
ciedlałby wszystkie szczegóły struktury pasmowej, W świetle wyników 
przedstawionych w tym rozdziale można sądzić,że za anomalny rozkład 
energetyczny ze ściany (100) wolframu w dużej mierze odpowiedzialna 


jest hybrydyzacja pasm 4 ;. Zarówno teoria, jak i doswiadczenie wske 
zują, że największe wzmożenie emisji umiejscowione jest w otoczeniu 
lokalnego maksimum dolnego złysrycyzowanego pasma, Szerokość 

tego pasma, uzyskana doświadczalnie, może różnić się od obliczonej 
teoretycznie. Spowodowane to jest wpływem pewnych efektów, związa- 


nych z anizotropowością struktury energetycznej wolframu, 


FIŚMIENNICTWO 


1 Slater WG, Koster G. R: Phys. Rev, 94, 1498 
(1954). 

2Mueller F, M.: Phys. Rev. 153, 659 (1967), 

3 Pltat M,Taranko R, Wysokiński K. L: 
Phys. Stat, Sol, (b) 71, K169 (1975), 


4 Taranko R; Phys. Stat. Sol. (b) 82, K5 (1977). 


PESDME 


3 paóoTe cOcyxaseTCA TpcÓleMa SBTOSIEKTPOHHO" MICHU 
K3 METSJAJIOB, KOTODHX 30HH8A CTDYKTYPS OTMCHHA C HOMOLIBO ABYX 
PHÓPAZUSYDANXCA 3CH. TolyueHWe SHEpreTHYecKUe pasTpeheN1eHHA 
3LETVUDOB3HHHX SNEKTPOHOB OTDAXADT BCE XSDBKTEDHPE TOUKU BOHHOŻ 
CTDYKTYDH 3MI:Tep3. 
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SU M MARY 


The problem of tno ficld emission of elec.rons from metals was 
investigated. The emitter was described by the broac anc narrow 
energy band, These bands were obtained in th: sinusoida! modet 
crystal potential. The total energy distribution of the emitted electrons 
shows all the interesting points of the band structure of the emitter. 


Złożono w Redakcji 10 VII 1978 roku. 


